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We show that the renormalization factor relating the renormalization group invariant quark masses to the bare 
quark masses computed in lattice QCD can be determined non-perturbatively. The calculation is based on an 
extension of a finite-size technique previously employed to compute the running coupling in quenched QCD. As 
a by-product we obtain the A-parameter in this theory with completely controlled errors. 



1. INTRODUCTION 

Calculations of the quark masses in lattice 
QCD are in principle straightforward, but there 
are several sources of systematic errors which 
must be carefully studied (see ref. |l| for a re- 
view of the status of these calculations and an 
up-to-date list of references). One of the uncer- 
tainties arises from the renormalization constant 
needed to convert from the lattice normalizations 
to the MS scheme of dimensional regularization. 
Usually one relies on bare perturbation theory (or 
some modified form thereof) to evaluate this fac- 
tor. Since only the one-loop term of the expan- 
sion is known, and since the gauge coupling is 
not small on the accessible lattices, the associ- 
ated error is, however, difficult to estimate. While 
the extension of the perturbation expansion to 
the next order may be helpful at this point, it is 
quite clear that a non-perturbative determination 
of the renormalization factor will be required to 
remove all doubts on the reliability of the quark 
mass calculations in lattice QCD. 

The calculation of the quark mass renormal- 
ization factor discussed in this talk is based on 



a recursive finite-size technique which allows one 
to compute the scale evolution of the renormal- 
ized parameters and fields from low to very high 
energies. An uncontrolled application of pertur- 
bation theory can thus be avoided. The general 
strategy of the computation has already been de- 
scribed in ref. ^ and it is our aim here to report 
on the progress that has been made along these 
lines in the case of the quark mass renormaliza- 
tion (for an alternative approach to the problem 
see refs. §1). 

2. PCAC RELATION 

Quark mass ratios can be accurately estimated 
using chiral perturbation theory (see ref. |^] for 
a recent discussion). To determine the absolute 
values of the quark masses in any given renormal- 
ization scheme it is hence sufficient to compute a 
particular linear combination of them such as the 
sum of the up quark mass and the strange quark 
mass. A possible starting point then is the PCAC 
relation 



9 ai («7m75s)r = (m u + m s )(u-y 5 s)n 



(1) 
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between the renormalized AS* = 1 axial current 
and the associated renormalized pseudo-scalar 
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density. 

In 0(a) improved lattice QCD the axial current 
and density are given by || 

(m7m75s)r = Z a (1 + b A am q ) (w7 M 7 5 s)iat, (2) 

(m75s)r = Z F (1 + b P am q ) (w7 5 s)i a t, (3) 

where (. . .)i at denotes the unrenormalizcd im- 
proved fields and m q is the average of the bare u 
and s quark masses (with the additive renormal- 
ization taken into account). Eq. (1) then holds 
up to terms of order a 2 . In the present context 
the factors 1 + bxam q amount to corrections of a 
few percent at most. The one-loop expression ]?j 

b A - bp = -0.001 x gl + ... (4) 

and a recent non-perturbative calculation of this 
difference || moreover suggest that these factors 
nearly cancel in eq. (1). At the present level of 
accuracy it thus appears safe to drop them. 

So if we determine (m u -I- m s )i a t through the 
vacuum-to-kaon matrix element of the unrenor- 
malizcd PCAC relation, 

9 /1 (tt7/ 1 75«)i a t = {m u + m s )i at (u75s)iat, (5) 
it follows that 

m u + fn s = (m u + m s )i at Z A /Z P . (6) 

The calculation of (m u + m s )i a t is standard by 
now |lj and will not be discussed here. As for 
the renormalization factors we note that Z A (g Q ) 
has been computed non-perturbatively || to a 
precision of about 1%, using a variant of the chiral 
Ward identity method of refs. fuJ-Q. 

The remaining unknown factor in eq. (6) thus 
is the renormalization constant Zp(go, afi). This 
factor is much harder to determine than Z A (go), 
because it is a function of two variables, the gauge 
coupling and the normalization scale (i of the cho- 
sen renormalization scheme. 

3. SCALE DEPENDENCE 

In the continuum theory the scale evolution of 
the running coupling and quark masses is gov- 
erned by the renormalization group equations 

Mt^ = 0(g), = T {g) m - (7) 
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Figure 1. Scale evolution of the running quark 
masses in quenched QCD 

For the case of dimensional regularization with 
minimal subtraction the 0- and r-functions have 
been calculated up to 4-loop order of perturba- 
tion theory ((l5|-|(l^| and references cited there). 
Using these results the evolution equations can 
be integrated starting at some large /i where the 
value of the running coupling is known. If we 
define the renormalization group invariant quark 
masses M through 

M = lim m(2b g 2 y do/2bo , (8) 

where bo and do are the one-loop coefficients of 
the (3- and the r-function respectively, the curves 
shown in fig. 1 are thus obtained. 

It should be evident at this point that the 
knowledge of the running quark masses at high 
energies is equivalent to providing the values of 
the corresponding renormalization group invari- 
ant quark masses. In other words, instead of 
m u + fn s we may just as well compute M u + M s . 
It is easy to show that the renormalization group 
invariant masses are non-perturbatively defined 
and that they do not depend on the chosen renor- 
malization scheme. We can hence compute them 
in any scheme that we like such as the one de- 
scribed in the following paragraphs. 

4. SF SCHEME 

The Schrodinger functional (SF) scheme is a 
particular finite-volume renormalization scheme 
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A k {x) = C k 



space 

(LxLxL box with periodic b.c.) 

Figure 2. Space-time manifold used to set up the 
SF scheme 

which has previously been used to compute the 
running coupling in the pure SU(3) gauge theory 
(l8[ . Its application in the present context has 
already been discussed in ref. [^) and in this pre- 
liminary report we only describe some of the key 
features of the scheme. 

The basic idea is to consider QCD in a fi- 
nite space-time volume of size L in all directions. 
Renormalization conditions are then specified at 
scale fi = 1/L and vanishing quark masses. To 
be able to perform numerical simulations at zero 
quark masses the boundary conditions are cho- 
sen in such a way that a frequency gap of or- 
der 1/L is induced on the quark and gluon fields. 
This can be achieved by assuming the space-time 
manifold to be as shown in fig. 2 and by impos- 
ing Dirichlet boundary conditions at time xo = 
and xo = L. More precisely, the spatial compo- 
nents of the gauge field are required to be equal 
to some prescribed constant abelian fields C and 
C . The response of the system to a change of 
these boundary values may be used to define a 
running coupling g 2 . We do not give any further 
details here, because the definition of g 2 that we 
have used is exactly the same as in ref. jTsfl . 

The renormalized quark masses in the SF 
scheme are defined through the PCAC relation 
and the renormalization condition for the pseudo- 
scalar density (cf. sect. 2). The precise form of the 
latter is of only practical importance since all def- 
initions lead to the same values of the renormal- 
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Figure 3. Quark diagrams contributing to fp(z) 
(diagram a) and fi (diagram b). The horizon- 
tal lines indicate the boundaries of space-time at 
xo = and xq = L. The pseudo-scalar density is 
inserted at the point z 

ization group invariant quark masses. A simple 
choice for the renormalization constant is 



Zp(g ,L/a) = cyffi/ f P (z)\ Zo=L/2 



(9) 



where fp{z) and f\ are unrenormalized correla- 
tion functions involving the pseudo-scalar density 
and the quark fields at the boundaries of space- 
time (c is to be chosen such that Z-p = 1 at 
go = 0). The corresponding quark diagrams are 
shown in fig. 3. Note that the renormalization 
factors associated with the boundary quark fields 
cancel in the ratio (9). In particular, the renor- 
malized coupling and quark masses defined in this 
way satisfy the usual renormalization group equa- 
tions [eq. (7)] with /i = 1/L and the appropriate 
[3- and r- functions. 

5. RENORMALIZATION GROUP 

In the SF scheme a change of the normaliza- 
tion scale amounts to a change of the lattice size 
at fixed bare parameters. By considering pairs of 
lattices with sizes L and 2L, we can thus study 
the evolution of the running coupling and quark 
masses under changes of the normalization scale 
by factors of 2. The important point to note is 
that the box size L can be as small as we like 
in physical units. The only restriction is that L 
and the low-energy physical scales in the theory 
(such as the radius r |il| ) should be significantly 
greater than the lattice spacing to avoid large cut- 
off effects. Such studies can be carried out us- 
ing numerical simulations and the scale evolution 
of the renormalized parameters is thus obtained 
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Figure 4. Simulation results for the step scaling 
function a(g 2 ) (quenched QCD) 

non-perturbatively. 

Close to the continuum limit the evolution from 
size L to 2L is described by the "step scaling func- 
tions" a and erp through 



f{2L) = a(f(L)), 
Z P (2L)=a P (g 2 {L))Zp(L). 



(10) 
(11) 



It is our experience that the lattice corrections 
to these equations are small and can be extrapo- 
lated away by repeating the computation of the 
step scaling functions for various lattice spacings 
at fixed g 2 . The step scaling function associ- 
ated with the coupling has first been calculated 
in ref. |Q and further data have since then been 
added so that this function is now known very 
accurately over a large range of couplings (see 
fig. 4). As for the other function, ap, the available 
data (fig. 5) are already quite usable and further 
runs will be made to fill the gap at large couplings 
and to achieve a better precision. All these results 
refer to the continuum limit, i.e. they involve an 
extrapolation to a = as indicated above. 

Once the step scaling functions have been de- 
termined, the sequence of couplings and quark 
masses 
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Figure 5. Preliminary results for the step scaling 
function ap(g 2 ) (quenched QCD) 



is obtained recursively using 

u k+ i = cr(u fc ), m k+1 = m k /ap(u k ). 



(13) 



u k 



f(2 k L), 



n^. = m(2 k L), 



(12) 



Since a(u) is monotonic the recursion can also be 
solved backwards, i.e. we can move up and down 
the energy scale as we wish. 

A technical problem here is that the step scal- 
ing functions are only known at certain values of 
the coupling and only to a finite numerical preci- 
sion. This difficulty can be resolved by fitting the 
data with a polynomial, as shown in the figures, 
and using the fit functions in the recursion (13). 
We have verified that the systematic error which 
is incurred by this procedure is neglible compared 
to the statistical errors provided one stays in the 
range of couplings covered by the data. 

6. RESULTS 

The largest value of g 2 which can be reached 
with the available data for the step scaling func- 
tions is 3.48. This corresponds to a certain box 
size -L max , which is large enough that contact can 
be made with the physical low-energy scales in the 
theory. Following ref. jl8| , and using some recent 
results of the UKQCD collaboration |f2fj| J2~l| , one 
finds L m ax/^o = 0.680(26). This can be taken as 
initial condition for the non-perturbative evolu- 
tion of the coupling and the recursion (13) then 
yields the data points shown in fig. 6. For com- 
parison we also plot the curves that one obtains 
by integrating the 2- and 3-loop evolution equa- 
tion starting at the right-most data point. Evi- 
dently the exact evolution of the SF coupling is 
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Figure 6. Evolution of the running coupling 
a(fi) = <7 2 /47r, fj, = 1/L, in quenched QCD, tak- 
ing ro = 0.5 fm to convert to physical units 

accurately matched by perturbation theory down 
to surprisingly low energies. 

We now digress a little bit and discuss how the 
A-parameter can be extracted from the calculated 
values of the running coupling. The idea is to 
evaluate the exact expression 



A = ^b g 2 )- bl/2b °e- 1 ^ 2 



x exp 



dg 



1 1 hi 



(14) 



at the high-energy data points of fig. 6. In this 
range of couplings the integral in the last fac- 
tor may be reliably calculated using the pertur- 
bation expansion of the /3-function, which (in 
the SF scheme) is now known to 3-loop order 
22 1 . Higher-order corrections are negligible at 
this point, and after converting to the MS scheme 
the result 



Ai21 = 251±21MeV 

MS 



(15) 



is obtained, where ro = 0.5 fm has been used to 
set the scale (the index (o) reminds us that this 
number is for quenched QCD). 

Now that the evolution of the coupling is un- 
der control, it is straightforward to obtain the 
data points in fig. 7 by solving the quark mass 
part of the recursion (13). Initially the recursion 



Figure 7. Evolution of the running quark mass 
given in units of the renormalization group in- 
variant mass M [eq. (8)] 

yields the ratio m(L)/m(2L max ). As in the case 
of the A parameter, the evolution may then be 
continued to infinite energy using the 2-loop ex- 
pression for the r-function in the SF scheme. In 
this way one is able to extract the renormaliza- 
tion group invariant mass M and hence the ratio 
m(L)/M plotted in fig. 7. Note that here too 
the scale evolution is accurately reproduced by 
perturbation theory down to the lowest energies. 
In particular, there is little doubt that the use 
of the perturbative evolution at high energies is 
safe. We would like to emphasize, however, that 
the absence of large corrections to the perturba- 
tive evolution of the running quark mass depends 
on the chosen scheme and should not be taken as 
a general feature of the theory. 

At the lowest-energy data point in fig. 7 we 
have 



L = 2L r , 



M/m = 1.18(2). 



(16) 



Recalling our discussion in sects. 2 and 3, it fol- 
lows from this that the total renormalization fac- 
tor, relating the bare current quark mass appear- 
ing in eq. (5) to the renormalization group invari- 
ant mass, is given by 

M/mut = 1.18(2) xZ A (. go )/M.9o,2L m ax/a). (17) 

The factors Za and Zp on the right-hand side 
of this equation are listed in table 1 for two val- 
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Table 1 

Evaluation of eq. (17) 



a 


2^ max/ & 




Z A 


M/mut 


6.0 


9.03(3) 


0.490(2) 


0.791(9) 


1.90(4) 


6.2 


11.63(2) 


0.500(2) 


0.807(8) 


1.90(4) 



ues of the bare coupling (3 = 6/c/q. Taking the 
product then yields the desired renormalization 
factor. These numbers are still preliminary, but 
they show what can be expected to come out once 
our study has been completed. 

We finally mention that the one-loop formula 

M/miat = (26 05 <;)-4/ii {i _ o.l2 x g 2 + . . .} (18) 

evaluates to 1.81 and 1.84 at = 6.0 and = 6.2 
respectively. This compares quite well with the 
non-perturbative result quoted in table 1. Note, 
however, that significantly lower values would be 
obtained if the expansion was written in terms of 
Parisi's boosted bare coupling. 

7. CONCLUSIONS 

In this work we have demonstrated that scale- 
dependent renormalization factors can be calcu- 
lated non-perturbatively without compromising 
approximations. The methods that we have em- 
ployed are completely general and are hence ex- 
pected to be useful in other contexts as well. 

Quark mass values are usually quoted in the 
MS scheme of dimensional regularization at a nor- 
malization mass /« = 1 GeV or [i — 2 GeV. 
Once a non-perturbative solution of the theory 
becomes possible, this convention is not entirely 
satisfactory, because the MS scheme is only mean- 
ingful to any finite order of perturbation theory. 
The renormalization group invariant masses, on 
the other hand, are non-perturbatively defined 
and scheme- independent. These quark masses are 
hence more quotable than the MS masses and we 
would like to recommend their use in future stud- 
ies. 

Previous computations of the A-parameter in 
LQCD have shown that this is a rather elusive 
quantity. The basic problem is that it refers 
to the high-energy limit of the continuum the- 
ory which is not readily accessible on the lat- 



tice. Applying our recursive method we have now 
been able to overcome this difficulty and to ob- 
tain the A-parameter with completely controlled 
errors (apart from quenching which remains to be 
the major limitation of present-day LQCD). 

This work is part of the ALPHA collaboration 
research programme. We thank DESY for allo- 
cating computer time to this project. 
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